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Abstract: We investigate the Sakai-Sugimoto model at nonzero baryon chemical potential 
in a background magnetic field both in the confined phase and in the deconfined phase with 
restored chiral symmetry. In this case the 8-brane Chern-Simons term becomes important. 
In the confined phase it generates a gradient of the pseudo-scalar "pion", which carries a 
non-vanishing baryon charge. Above a critical value of the chemical potential there is a 
second order phase transition to a mixed phase which includes also ordinary baryonic matter. 
However, at fixed baryon charge density the matter is purely "pion" -gradient above a critical 
magnetic field. In the deconfined chiral-symmetric phase at nonzero chemical potential the 
magnetic field induces an axial current. We also compute the magnetization of the baryonic 
matter and find that it is paramagnetic in all three phases but with nonlinear behavior at 
large magnetic field. 
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1. Introduction 

The behavior of QCD under external conditions is an interesting and physically relevant 
problem. For example, at high temperature the ground state is believed to be a deconfined 
quark-gluon plasma, and at high density it is believed to be a color-superconductor. The for- 
mer is relevant for the understanding of the physics at RHIC, and the latter may be relevant 
for the physics of dense stellar objects such as neutron stars. Background electromagnetic 
fields provide another kind of external condition. Their effect on the QCD ground state, in 
particular in the charged flavor sector, is possibly relevant for magnetars, which are neutron 
stars with very large magnetic fields. However, it is often the case that in the physically 
relevant regimes QCD is strongly coupled, and we cannot reliably use perturbation theory. 1 
Lattice gauge theory has been very successful for studying equilibrium properties of QCD at 
nonzero temperature, but it doesn't do so well at nonzero chemical potential, and it is not 
equipped to handle real-time transport properties, such as those expected in a background 
electric field. Another possibility which is non-perturbative is to use some truncation of 
the Schwinger-Dyson equations to study QCD. An alternative approach to strongly-coupled 

1 This appears to be the case at RHIC, for example. 
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gauge theories has emerged in recent years from string theory via the AdS/CFT correspon- 
dence and its generalization to gauge/gravity holographic duality. The holographic approach 
is well-poised to address questions related to external conditions since these translate simply 
to boundary conditions on internal fields in the bulk. For example, the temperature corre- 
sponds to the size of the Euclidean time dimension at the boundary, and a chemical potential 
associated with a conserved current corresponds to the boundary value of the bulk gauge field 
dual to the current. Determining how an external condition affects the ground state of the 
field theory is then just a matter of solving the bulk equations of motion with the appropriate 
boundary conditions. 

Due to asymptotic freedom, the holographic dual of QCD cannot be a (super)gravity 
theory alone and must include all the closed string excitations. Neverthless, it may be useful 
to study gravitational models as holographic models of low-energy properties of large N c QCD, 
in the hope of eventually embedding them in string theory with the correct UV properties. 
This has been the approach of the so-called "bottom-up" models. The "top-down" approach, 
in contrast, is to consider the full string theory in a background where the low- lying excitations 
resemble those of QCD and in which one can consistently study the supergravity limit. This 
does not give QCD, but these kinds of models share many of its low-energy properties. It is 
worth emphasizing again that we are studying these models in the hope that they resemble 
QCD, but the analysis is of these models only. The closest so far to QCD is the Sakai- 
Sugimoto model Q. This model consists of N c 4-branes wrapping a circle with anti-periodic 
boundary conditions for the fermions, Nf 8-branes at a point on the circle, and Nf anti-8- 
branes at another point on the circle. The low-lying open string excitations are precisely 
those of SU (N c ) Yang- Mills theory with Nf flavors of massless quarks. The holographic limit 
corresponds to N c — > oo, and Nf is kept finite so that the 8-branes are treated as probes in 
the near-horizon background of the 4-branes. In this limit the background is capped off in 
the IR, which corresponds to confinement in the gauge theory, and the 8-branes and anti-8- 
branes connect into U-shaped 8-branes, which corresponds to chiral symmetry breaking in 
the gauge theory. This model has also been studied in various external conditions, including 
nonzero temperature [§], nonzero baryon chemical potential H ||], and background electric 
and magnetic fields ||, |(| f7), in which it exhibits many properties that are expected of QCD. 

Background magnetic fields are particularly interesting in that they may be physically 
relevant in neutron stars, where they can reach values of about 10 15 Gauss. On the theory 
side, background magnetic fields have some interesting effects on the QCD ground state. One 
effect is the catalysis of chiral symmetry breaking by a strong magnetic field H|. The basic 
mechanism for this is that in a strong magnetic field all the quarks sit in the lowest Landau 
level, and the dynamics are effectively 1+1 dimensional. The effect of the magnetic field on 
the quark condensate was studied in ||]. The effect of a background magnetic field in the 
Sakai-Sugimoto model was studied in 0, where the catalysis of chiral symmetry breaking 
was demonstrated explicitly. In particular it was shown that the critical temperature for 
the restoration of chiral symmetry increases with the magnetic field and approaches a finite 
temperature at infinite magnetic field. 
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In this paper we will be interested in the effects of a background magnetic field at nonzero 
baryon chemical potential. This question was recently studied in the the low-energy effective 



field theory |i0| , 11, p^| , where it was shown that the triangle anomaly leads to interesting 
effects. 2 In the deconfined chiral-symmetric phase the combination of a magnetic field B and 
a nonzero baryon chemical potential hb leads to a non-zero axial current density [|TTf 



3a = y^HbB. (1.1) 

This current is generated purely by fermionic zero modes and is therefore topological in 
nature. The result is therefore also exact. In the confined phase the anomaly leads to a 
non-trivial pion gradient and an associated baryon charge density fljjfl, 



^=^bB, i=^B.V,". (1.2) 

We will show that similar effects occur in the one-flavor Sakai-Sugimoto model. In this 
model the anomaly is encoded in the five-dimensional Chern-Simons term of the 8-brane 
action. The model does not include a true electromagnetic field, but we can mimic the effect of 
a non-dynamical (background) electromagnetic field using the non-normalizable mode of the 
8-brane world volume gauge field. This field is actually dual to the baryon current in the gauge 
theory but is equal in the one-flavor case to the electric current. We will therefore use the same 
bulk gauge field, but different components, to describe both the baryon chemical potential 
and the background magnetic field. We will show that these source a third, normalizable 
component of the gauge field via the Chern-Simons term. In the low-temperature confining 
background this field has a nonzero boundary value, which is interpreted as the gradient 
of the U{\)a pseudo-scalar meson, i.e. the rj'. This also leads to a baryon number charge 
density. For small magnetic fields, our result agrees with ( |l.2[) adapted to the U(1)a sector. 
Furthermore, we will show that there is a phase transition at a critical value of the chemical 
potential to a mixed phase of ordinary baryonic matter and pseudo-scalar gradient matter. 
In the mixed phase the relative proportion of ordinary baryonic matter at fixed chemical 
potential decreases with the magnetic field. In the high-temperature deconfining background, 
in the restored chiral-symmetry phase, the induced gauge field has a vanishing boundary 
value, and the leading asymptotic behavior corresponds to an axial current density, which 
agrees with (|Q]), 

We will make two simplifying assumptions about the 8-brane embedding, which do not 
affect our results qualitatively. First, we will consider only the one-flavor case, in which the 
8-brane worldvolume gauge field is abelian. This will allow us to use the full DBI action for 
the 8-brane. Second, we will consider only the antipodal 8-brane embedding. This means 
that in the low-temperature confining background the tip of the 8-brane coincides with the 
tip of the space. This will also simplify the analysis with sources since the embedding will 
remain smooth. This should not affect the qualitative results since this embbeding is smoothly 



2 For other interesting anomaly- induced effects at finite density and magnetic field see jl3 
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connected to the non-antipodal embedding. 3 In the high-temperature deconfining background 
it implies that the preferred embedding is the disconnected 8-brane-anti-8-brane configuration, 
in which chiral symmetry is restored. Finally, to avoid clutter we will work mostly with 
dimensionless quantities by absorbing appropriate factors of R and a'. We will denote such 
quantities by lower case letters. For example the dimensionless coordinates are u = U/R 
and x^ = Xu/R- A translation table between the dimensionless and physical versions of the 
relevant quantities is provided in the appendix, but we will also give the translation whenever 
a new quantity is introduced. 

The rest of the paper is organized as follows: In section 2 we review the relevant features of 
the Sakai-Sugimoto model at nonzero baryon chemical potential. In section 3 we analyze the 
magnetic properties of the confined phase, including the pseudo-scalar gradient and baryon 
charge density, as well as the magnetization. In section 4 we study the magnetic properties 
of the deconfined phase. Section 5 contains our conclusions. 

Note: While our paper was being completed the paper [14], with which there is some overlap, 
came to our attention. 

2. Review of finite density HQCD 
2.1 Confined phase 

The Euclidean background dual to the confining phase is given by 

d.S 2 = U" 2 {{dx$) 2 + dx 2 + f{u)dx\) + U- 3 / 2 (j^ + U 2 dn 2 ^ 

-* - g s u 3/4 , F 4 = 3ir(a') 3/2 N c dn A , (2.1) 



e 



where X4 is a compact coordinate with periodicity 2ttR a , and 

/(.) = !-%, u KK = § t . (2.2) 

The curvature radius of the space is given by 

R = {*g,N c fl*>/d , (2.3) 
and this is related to the four-dimensional 't Hooft coupling 



A = * — . (2.4) 

The antipodal embedding of the 8-brane in this background has a U shape that satisfies 
x'^iu) = 0, with the tip at ukk- Other than the embedding scalar field X4, the 8-brane 
worldvolume theory contains fermions and a gauge field. We will ignore the fermions. The 



3 in fact there is a scaling argument connecting all of these embbedings, see the second reference in 
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gauge field has, in general, both a vector and an axial part depending on the parity under 
exchanging the two halves of the embedding, 



a M (x», u) = a v M {x», u) + a^x", u) . (2.5) 

The physical gauge field is Am = aM-R/(27ra'). We will consider only fields that are uniform 
on the S 4 , so the index // runs over — 3. There is a discrete spectrum of normalizable 
radial modes corresponding to various low-spin mesons. In particular the zero mode of 
is identified with the massless pseudo-scalar corresponding to the Goldstone boson of the 
broken chiral symmetry. For a single flavor this is the rf . However there is some freedom 
in identifying the mesons due to the gauge symmetry. A particularly nice gauge choice, that 
preserves the four-dimensional Lorentz symmetry, is a u = In this gauge the pseudo- 
scalar reappears in the zero mode of a^, 

a^(x^, u) = d^ip(x^)ipo(u) + higher modes , (2.6) 

where 




ipo(u) = — arctan t — ^ 1 . (2-7) 



The physical pseudo-scalar field is rj'(X) = f n p(x)R 2 /(2ira'), where f n is given by 5 



3/2 



2 N c u Kjx 

f * = ~8^r- (2 ' 8) 

Note that the axial zero mode has a normalizable field strength. By contrast, the zero mode 
of the vector part of the gauge field a^{x^) is n-independent and therefore non-normalizable. 
In general, this corresponds to a source for the vector (baryon) current in the boundary gauge 
theory. In particular, the asymptotic value of the ^-independent part of Oq is identified with 
the baryon chemical potential, 

a o ( u ~* °°) = M • (2-9) 

In our convention the baryon charge of a quark is 1, rather than 1/N C . 

For a static and uniform baryon charge distribution the (Euclidean) 8-brane DBI action 
per unit 4-volume of spactime is given by 



S DB i=M [°° duu 5 / 2 x U--(aV>( U )y, (2.10) 

Ju KK y f(u) 



4 The anomalous mass of the rj is suppressed at large N, 

5 The pion decay constant was determined in terms of the parameters of the model by comparing the 
non-abelian 8-brane Yang-Mills action with the standard non-linear sigma model 
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where the overall normalization is given by 

AT E>2 

The factor of 2 corresponds to the two halves of the embedding. The resulting equation of 
motion for the gauge field is given by 



d_ 

du 



yi-/(«)(c#'(«)) 2 . 



0. (2.12) 



Integrating once gives 



<#'(") = 1 * (2.13) 

\/f{u) Vr + cr 



where d is the constant of integration. The asymptotic solution at large u is then 

(2 - 14) 

Since the action (per unit 4-volume) of the solution defines the grand potential (per unit 

3- volume) of the gauge theory at a fixed /i, we identify the constant d as the baryon charge 
density. The physical chemical potential is /ig = /J,R/(2ira'), and the physical baryon charge 
density is D = d(2ira'N / R). 

In the absence of sources the only solution is a constant 

al(u) = n- (2.15) 

In this case the gauge field is pure gauge, and therefore the physics does not depend on 
the value of [/,. A second solution becomes possible when one includes sources at the tip 
corresponding to 4-branes wrapped on the S 4 . These 4-branes are precisely the baryons of the 
model. A single baryon carries N c units of baryon charge. Assuming a uniform distribution of 

4- branes with positive number density 77,4, and assuming that the 4-branes are well-separated 
in space so that we can ignore interactions between them, the source action per unit 4-volume 
is given by 

<Sd4 = AA ^77,4777,4 — n^Nc j dua^ (u)5(u — ukk)^J ■ (2-16) 

The first term comes from the 4-brane DBI action, where 777.4 is the mass of a wrapped 4-brane, 
i.e. a baryon, located at u = ukk, 

m-4 = 7^N C UKK ■ {2.17) 

The physical 4-brane mass and density are given by M4 = ra^Rj (2ira!) and ./V4 = n^ftna'Af/R). 
The second term in the source action comes from the N c strings that connect each 4-brane 
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to the 8-brane (or equivalently from the 8-brane CS term, if we describe the 4-branes as in- 
stantons in the 8-branes p5| ). This relates the baryon charge density to the 4-brane number 
density as 

d = N c n A . (2.18) 

We can then determine this number for the solution by extremizing the action with respect 
to n.4. This gives a condition on the gauge field at the tip, 6 

<%{u KK ) = j± (2.19) 

which implies that the solution exists only for /x > ra^jN c ? This has the obvious phenomeno- 
logical interpretation that, at low temperature, baryons can only appear when the chemical 
potential is high enough to produce them. Furthermore, it is easy to see from the form of 
the action that this "nuclear matter" solution dominates over the vacuum solution ( [2.15 ). 
In other words, nuclear matter forms as soon as it can, and there is a phase transition at 
\x c = m^jN c . The relation between the baryon charge density and chemical potential is 
obtained by integrating (2.13), 

f°° du d . 

Ju KK v/W v u +« 
Near the critical point we get a linear relation 



3u 



3/2 
KK 



dn—™^-^), (2.21) 



7T 



which implies that the phase transition is marginally second-order. This is different from 
the expected first-order transition in QCD. However the result is reasonable since we have 
ignored baryon interactions. 

2.2 Deconfined phase 

The Euclidean background for the deconfined phase is given by 

ds 2 = J {f(u)(dxf) 2 + (fx 2 + dx\) + «-f (J^ + u 2 d^ , (2.22) 
with the same dilaton and RR field as before, and with 

/(«) = 1 - 4 > (2-23) 

where ut is related to the temperature, i.e. the inverse periodicity of the Euclidean time Xq, 
as ut = (4tt/3) 2 R 2 T 2 . For T > l/(2vri? 4 ) this phase dominates over the confined phase and 



6 This can also be seen by requiring a consistent interpretation of the themodynamic potentials 0| 
7 For anti-four-branes this is ~m±jN c . 
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the the theory deconfines. In this phase the 8-branes in general have two possible embeddings: 
a U-shaped embedding similar to the one in the confined phase, and a parallel 8-brane-anti- 
8-brane embedding M. We will consider only the parallel embedding, which is the dominant 
phase at all temperatures in the antipodal case. 8 In this embedding there are two independent 
gauge fields and a^, associated with the 8-brane and anti-8-brane, respectively. As before 
we work in the gauge a u = a u = 0. 

There are a number differences from the confined phase. First, the spectrum of normaliz- 



able solutions is not discrete, so there is no particle (meson) interpretation [16|. Second, both 
zero modes are non-normalizable and therefore correspond to two sets of parameters in the 
gauge theory. We will set the axial parameters to zero. Our boundary conditions at infinity 
are therefore 

ao(oo) = ao(oo) = \x . (2.24) 

Another important difference is that for regularity at the horizon we must impose the bound- 
ary conditions 

a (u T ) = a (u T ) = 0. (2.25) 

The two gauge fields are therefore equal, and the total action for the 8-brane and anti-8-brane 
is given by 

S DB i=M [ duu 5/2 ^l-{a' (u)) 2 , (2.26) 



where the normalization is the same as in (2.11), with the factor of 2 accounting for the two 
branes. The solution now satisfies 9 

f°° d 

H = a (oo) = / du . (2.27) 

J ux v u ~\~ d 

In this phase matter is made up of deconfined quarks and begins to form immediately at 
nonzero /j,. For small // the density is given by 

dra — J-/i~T 3 /i- (2.28) 

3. Magnetic properties of the confined phase 

To mimic the effect of a background magnetic field we turn on a background value for the 
zero mode of a spatial component of the vector gauge field, 

a^(x2,u) = hx 2 . (3.1) 



s Below a certain value of the asymptotic 8-brane-anti-8-brane separation there is a range of temperatures 
for which U-shaped embedding dominates, and the theory realizes an interesting intermediate phase of decon- 
finement and chiral symmetry breaking. We will not consider this phase here. 

9 This can be expressed in terms of a hypergeometric function ]1 7|] . 
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The physical magnetic field is H = h/{2ira.'). Since %(u) 7^ 0, the five-dimensional CS term, 
which comes from the 8-brane CS coupling to F4, will source the axial field af(u). As we 
saw in the previous section, the boundary value of this field corresponds to the (constant) 
gradient of the pseudo-scalar (in the x l direction in this case), 



Vtp. 



(3.2) 



This corresponds to a field, rather than an external parameter, in the gauge theory, since the 
zero mode of the axial field is normalizable. Its value is therefore determined by extremizing 
the action. Furthermore, since is an axial field, it must vanish at the tip, 

o>i{u K k) 



0. 

The action with all the relevant fields is Sdbi + Scs, where 10 
dun^ 2 



Sdbi 
Scs 



00 




{aV>{u)Y + (af(u)) 



1 + 



h 2 



(3.3) 



(3.4) 



-M I du (8203 Oq (u)a 1 '(u) — <92«3 %' \ u ) a i ( u ) 

'U-KK 



Va V Ai , Va A V>\ 
a 3 d 2 a a x + a| d 2 a{ a v ) , 



The corresponding integrated equations of motion are given by 

Vu 5 + h 2 u 2 a%'(u) 



1 

7R 



K u)) 2 + (af(u)f 
Vu 5 + h 2 u 2 af(u) 



7R 



3haf(u) + N c n 4 



3/iag (n) + c , 



(3.5) 



(3.6) 



where the constant of integration in the equation has been identified with the density of 4- 
brane sources as before, and the constant of integration in the equation will be determined 
shortly. 

The action as written above is a bit problematic however. As is well known, in the 
presence of a boundary the CS action is not invariant under gauge transformations that do 
not vanish at the boundary, and requires the addition of a boundary action whose gauge 



transformation cancels that of the bulk action 18]. This will be true also in the infinite 
volume case if we allow for gauge transformations that do not vanish at infinity, which we do. 
In our case the gauge is partially fixed by a u = 0,2 = 0, and by the identification of a^(oo) 
and af(oo) with ju and V99, respectively. This leaves a residual gauge symmetry that depends 
only on X3, under which — ► + dsA(xs). Under this transformation the CS action above 
transforms by (restoring the integration over the x^, and integrating by parts symmetrically): 



5 Scs 



2 



dxodxidxs dsAfas) 





X2—+00 






J dua^a^' 




[ dx2d 2 a^a^ 




X2^ — OO •> 


u^oo 



(3.7) 



10 Our ansatz is that ao(u) and ai(w) depend only on u. However, we retain two terms in the CS action with 
derivatives with respect to X2, since they contribute to the equations of motion. 
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In particular, the first term will be non-vanishing on shell since the gauge transformation 
does not depend on x 2 - Subsequently we need to add a boundary action at the boundaries 
X2 — > ±00 and u — * 00 that has the form 



^bndy 



" 2 



f ,4 V a V A 
+ Id xa 3 o 2 a^ a^ 



(3.8) 



After integrating by parts the last two terms in Scs, the sum of the CS and boundary actions 
becomes 



/T3 1 
d A xdu -d 2 a%a[ a^ - -a\'aJ Q d 2 a^ 

which is manifestly invariant under the residual gauge transformation. 11 
The baryon charge and currents are defined by 



(3.9) 



J" Or) 



(3.10) 



5A /1 (x, u = 00) 

Where S eom is the value of the action on the equation of motion. This can be computed by 



6S, 



ESC ro 5C 
T7r-;bdiA + —5A 



6d { A 



5A 



(3.11) 



By integrating by parts and using the equation of motion we find 

5C 



J"(x) = lim 

u^oo \dd u A^{x) 

where the right hand side is evaluated on the equation of motion. 
We can now read off the (dimensionless) baryon charge density 

3 3 
d = N c n& + — ha-i (00) = N c n^ + —hX/p . 



(3.12) 



(3.13) 



The origin of the second term can be understood as an additional 4-brane charge inside the 
8-brane, which is due to the orthogonal worldvolume field strengths in the (x 2 , x 3 ) and (u, x 1 ) 
directions. We can likewise get the (dimensionless) axial current density from (|3.6|) , 



3 3 
3A = c + -ha% (00) = c+ -hfj, . 



(3.14) 



Recall however that we still need to extremize the action with respect to V93. This has the 
effect of setting ja = 0, and therefore c = —\h[L. 



11 Note that with this modification the four dimensional current densities can be expressed as integrals over 
the radial direction of the appropriatelly defined five dimensional current densities. The five dimensional 
current densities are defined in general only up to the addition of an arbitrary constant multiple of *(F A F). 
The restriction of the gauge symmetry above fixes this constant to be i. 
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We can simplify the equations of motion considerably as follows. First, dividing (|3.5| ) by 
( |3.6| ) gives an expression that can easily be integrated to give the relation 

\h{al{u)) 2 - hfiaUu) = \h(at{u)) 2 + N c n m af(u) + « , (3.15) 

where k is a constant that depends on the type of solution. As in the zero-magnetic field case, 
there are two types of solutions, with and without 4-brane sources. In the sourced case there 
is an additional condition on olq at the tip given by ( |2.19| ). Using the values at the boundary 
in the sourceless case, and at the tip in the sourced case, we get 



h{W(p) 2 sourceless case 



2 

( n - ) sourced case . 



(3.16) 



Next, define a new coordinate 



, 3hdu 

y = I 1 — ■ (3-17) 

K (l + ||) + (iV> 4 ) 2 - (fV) 2 - 6hK 

Using the relation (3.15), and some algebra, the equations of motion then reduce to 



al'{y) = af(y) + ^ (3.18) 

4'(y) = aX(y)-l, (3-19) 

where the derivative is with respect to y. Let us now analyze the two types of solutions. 
3.1 Pseudo-scalar gradient phase 

In the absence of sources = 0, and all the baryon charge density comes from the pseudo- 
scalar gradient 

d=-hVtp. (3.20) 
The solution to ( |3.18| ) and ( p. 19 ) is given in this case by 



a , , u sinh y 
4(V) = £ , 3.22 

2 coshyoc 

where yoo = y(u — > oo) can be determined numerically in terms of /j, and h from the integral 
equation 

y - f°° 3hdu (3 23) 

i% K Z 1 / 2 \/ + h 2 u 2 — h 2 ^ 2 sech 2 yoo 
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W ~ ^irrVBH (3.26) 



The pseudo-scalar gradient is then simply 

Vip = |tanh yoo . (3.24) 

The numerical results for Vip and d as functions of /j, and h are presented in fig. [l|. For /j<1, 
i.e. sub-string scale magnetic fields, the behavior is linear, and the pseudo-scalar gradient is 
approximately given by 

2u KK 

As /i increases the nonlinearity of the DBI action becomes apparent. 
In terms of the physical quantities we get 

N c 

for small magnetic fields. This agrees with the one-flavor version of the result (|1.2| ) from 
The relative factor of N c /2 is understood as follows. First, the baryon charge of a quark in 
Q is 1/N C so / u^ ere = N c ^ re . 12 Second, the CS term coupling the baryonic U(l)v gauge 
field to the C3-component of the SU(2) gauge field has a factor of 2 relative to the purely 
abelian CS term once the boundary term at spatial infinity is omitted. 

3.2 Mixed phase 

Above some value of the chemical potential [i c a solution with 4-brane sources (baryons) is 
also possible. In this phase both the baryons and the pseudo-scalar gradient contribute to 
the baryon charge density. We expect a phase transition to occur at this value of [i to the 
mixed phase. We will compute the critical value as a function of the magnetic field fJ. c (h), 
as well as the total baryon charge density d(h, //) and the fraction of the total baryon charge 
carried by baryons. 

Using the boundary conditions at the tip Oq (y = 0) = m^/Nc, af(y = 0) = 0, the 



solution to the equations of motion ( |3.18|) , (|3.19|) is now given by 



a o(v) = (jf ~ f J coshy + ^1 sinhy + | (3.27) 

4{V) = - f ) sinh y + ^ (cosh y - 1) . (3.28) 

The boundary conditions at infinity then determine the gradient and 4-brane density implic- 
itly in terms of \i and h, 

_ coshj/op - 1 m 4 ^ 2g ^ 

sinh t/oo N c 

\h[i + \h\n- ^MJ coshyoo 

N c n 4 = ^— C -L , (3.30) 

sinh j/oo 

12 Note that since f n ~ N c , the pseudoscalar gradient is suppressed at large N c at fixed f/g ere , as expected 
in an anomaly-mediated effect. 
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d Vn 




Figure 1: The baryon number density d and the pseudo-scalar gradient V</5 as functions 
of /i for fixed h = 1, and as functions of the magnetic field h for fixed fi = 0.2, all with 

UKK = 1. 



where is the solution to to the integral equation 

3hdu 




sinh j/ 



r»4 



+ 



^wf) ( cosh 



+ 1) 



(3.31) 



The critical value of the chemical potential corresponds to the point at which the actions 
of the and mixed phases are equal. But it also coincides, as it did in the zero magnetic 
field case, with the minimal value of the chemical potential to create a baryon. This can be 
determined numerically by setting 774 = in ( |3.30| ). The result is the phase diagram in the 
(/i, h) plane shown in fig. |2[ The critical chemical potential increases from its h = value 
m^/Nc to 27714/ N c as h — > 00. For a given h there is a marginally second-order phase transition 
from the V</? phase to the mixed phase at /j c (/i) that generalizes the ordinary nuclear matter 
transition at h = 0. We also see that for a fixed total baryon charge density the pseudo-scalar 
gradient phase dominates above a critical magnetic field. 

The total baryon charge density is given by 



, AT 3,_ 3h 
a = iV c n4 + -hvip = — 



772-4 



cosh Voo + l 
sinliT/oo 



(3.32) 



Figure y shows the total baryon charge density and the fraction of that charge carried by 
baryons, which are obtained by numerically computing ( 3.32 ) and ( 3.3C| ) . We see that the 
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Figure 2: Phase diagram in the (a) canonical and (b) grand canonical ensemble. 



relative proportion of baryons at fixed h increases with [i. In the limit of large [i the system 
is almost entirely baryonic nuclear matter. On the other hand at fixed (i the proportion of 
baryons decreases with h. For m^jN c < fj, < 1m^jN c the proportion of baryons vanishes at 
the critical magnetic field h c (n) shown in fig. [2], where the V(p phase takes over. 




Figure 3: The total baryon charge density d and the baryon fraction n±N c /d, as functions 
of fj, for fixed h = 1 and ukk = 1, and as functions of h for fixed /j, = 'Sm^/Nc and ukk = !• 



3.3 Magnetization 

The state described by either the pseudo-scalar gradient phase or the mixed phase responds 
to the external magnetic field by getting magnetized. The magnetization M can be defined 
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in either the grand canonical ensemble as 



M(jjt,h) 



dQ(fi, h) 



dh 



(3.33) 



where the grand potential Q(/j,,h) = S[ao(u), ai(u)]\ EOM , or in the canonical ensemble as 



M(d, h) 



dF{d, h) 



dh 



(3.34) 



where the free energy F = £1 + [id. The magnetic susceptibility describes the linear response 
of the system to small magnetic fields and is defined as 



X 



dM 



dh 



h=0 



(3.35) 



in either the canonical or grand canonical ensemble. 

We would like to focus here on the matter contribution to the magnetization and sus- 
ceptibility. The magnetic properties of the vacuum were studied in [||, ||]. We will therefore 
subtract from the quantities above the formally divergent contribution of the vacuum, which 
gives a finite result that represents the corresponding contribution of just the matter. The 
numerical results for the magnetizations in the pseudo-scalar gradient and mixed phases are 
presented in fig. ||. For small h the response is linear, but for h ~ 0(1) the non-linear effect 
of the DBI action becomes pronounced. 




h 




Figure 4: The magnetization M (in units of jV) as a function of h in (a) the pseudo-scalar 
gradient phase for fixed u = 0.2 and (b) the mixed phase for fixed /i = im^/Nc (red) and 
fixed d = 1 (dashed blue), all with ukk = 1- 



The magnetic susceptibilities can be computed in a similar way. In the phase the 
grand-canonical magnetic susceptibility can actually be determined analytically to be 

Ax = x(A0 " X(0) = • (3-36) 

4u KK 

The canonical susceptibility in this phase can be computed numerically, although we will 
not do this here. In the mixed phase, the magnetic susceptibility can only be computed 
numerically. Figure ||| shows the susceptibility in the mixed phase, in both the canonical 
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and grand canonical ensembles. Our results show that the matter is paramagnetic in both 
phases. Interestingly, the contribution to the susceptibility from the DBI term is negative, 
i.e. diamagnetic, in both phases. However the contribution of the CS term is always positive 
and larger. 




Figure 5: The magnetic susceptibility A% (divided by TV) of the mixed phase as a function 

of d and \x for ukk = !• 



4. Magnetic properties of the deconfined phase 

The 8-brane DBI and CS actions in the deconfined background are given by 

Sdbi = M j™ duu 5 / 2 ^ (l - (a(,H) 2 + (l + ^ (4.1) 

POO 

Scs = -M / {d 2 aXal{u)4'{u) - d 2 a% a%'(u)a?(u) - ag^a? + a^d 2 afa^) , 

where now / = 1 — (ii|,/ii 3 ), and where we have included both the 8-brane and anti-8-brane 
parts, with ao = «o and a\ = —a±. The boundary value of the axial field is now a parameter, 
rather than a field, in the gauge theory, which we set to zero, ai(oo) = 0. As before in the 
confined phase, we modify the Chern-Simons action by throwing away boundary terms of the 
form 

-d 2 (aza^d^ + ^d u (a 3 «9 2 a[ia ]) (4.2) 

to obtain the correct five-dimensional currents. 

The integrated equations of motion are then given by 



Vu 5 + h 2 u 2 a' (u) 

- 3hai(u) + d (4.3) 



2 



Vi-K(^ + /(u)K(n)) 

V^TWf(u)a[(u) =3hao{u)+jA _K^ (4 . 4) 
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where d is the baryon charge density and ja is the axial current density. Unlike in the confined 
phase, we do not get an additional condition by extremizing the action with respect to ai(oo). 
However there is an additional condition imposed by regularity at the horizon, ao(^T) = 0. 
Since f(ux) = as well, the consistency of the a\ equation of motion flh4| ) requires turning 
on a specific axial current density 



3A 



-hn. 



In terms of physical quantities, the axial current density is 

T Nc H 
A = 4^ H ^ B ' 



(4.5) 



(4.6) 



which agrees precisely with the result ( |1.1| ) of 1 1 ] , once we account for the different normal- 
izations (the relative factor of N c /2) as in the pseudo-scalar gradient case in the previous 
section. 

The coupled equations of motion can be solved numerically using a shooting algorithm. 
The results for /i(d, h, T) are shown in figure ^. We see that at a fixed h, [i grows linearly with 
d for small d. This is the expected behavior of free massless fermions in 1 + 1 dimensions, 
which one may think is a natural result for massless fermions in 3 + 1 dimensions in the 
background of a large magnetic field. It is interesting to note however that the linear region 
extends to d ~ h , which is beyond the linear region for free fermions that extends to d ~ h. 
The deviation is not surprising since the fermions are not free. We see also that, at a fixed d, 
fx decreases with h. This reflects the increase of the ground state degeneracy with h. 



Figure 6: The baryon chemical potential /i as a function of (a) d with fixed h = 1 and 
T = 0.3, (b) of h with fixed d = 1 and T = 0.3, and (c) of T with fixed h = 1 and d=l. 



The magnetic response of the deconfined phase is quite similar to that of the mixed 
confined phase. The magnetization and magnetic susceptibility are defined, as in the confined 
phase, by (|3D and fl335p , now with the susceptibility of the deconfined vacuum subtracted. 
The numerical results for Ax(d,T) and M(h) computed in the canonical ensemble for fixed d 
are shown in figures ^ and ||. In particular, the high-temperature behavior of the susceptibility 
at fixed density is x ~ ^/T 9 , which deviates from the Curie law x ~ 1/T. 

Finally, it is also interesting to note that the equations of motion ( |4.3| ) and (4.4) are 
almost symmetric under the interchange of axial and vector components. We considered a 
vector ao and an axial ax, with a chemical potential only for the vector charge /x = oo(oo), 
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Figure 7: The magnetic susceptibility Ax (divided by J\f) of the deconfined phase as a 
function of (a) d for fixed T = 0.3 and (b) as a function of T for fixed d = 1 . 




Figure 8: The magnetization M (divided by Af) as a function of the magnetic field h in 
the deconfined phase for fixed d = 1 and T = 0.3. 



leading to an axial current 2A- However, if we consider instead an axial ao and a vector a±, 
with an axial chemical potential /j,a = ao(oo), this would lead to a non-zero baryon number 
current jb = \h^A- 



5. Conclusions 

In this paper we have explored the properties of one-flavor holographic QCD at finite density 
in a background magnetic field. It turns out that this system has a rich phenomenology. In 
particular, in the confined phase turning on a magnetic field induces a gradient for the pseudo- 
scalar field. This gradient carries baryon charge, and at large enough magnetic fields it is the 
dominant phase. That is, if we start at zero field with some baryons, as we increase the field 
those baryons will start being replaced by a gradient of the rf field, eventually disappearing 
altogether. In the chiral-symmetric deconfined phase we found that the magnetic field induces 
an axial current whose value is independent of the temperature. The first property can be 
traced to the axial anomaly of fermions, and the second phenomenon can be traced (at weak 
coupling) to the existence of particular fermionic zero modes in a magnetic field background. 
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In the holographic dual both of these properties are induced by the Chern-Simon term on the 
8-brane, but the second is also due to the presence of a horizon in the spacetime geometry. 
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A. Dimensional translation table 



quantity 


dimensionless variable 


physical variable 


coordinates 




Xf, = Rxp, U = Ru 


gauge field 




A - -R-n 

A 4 ~~ 2na' U V 


magnetic field 


h 


H = ^h 


baryon chemical potential 






baryon charge density 


d 


D = 2 ™ R N d 


axial current density 


3A 


J A- R J A 


pseudo-scalar field 


f 




wrapped 4-brane mass 




M 4 = ^m 4 


wrapped 4-brane density 


Tl4 


AT 4 = 2 ™'^ n4 
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